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Abstract 

We compute the dimension-2 gluon pair condensate g 2 (AA) and the dimension-4 mixed 
quark-gluon condensate (qgA ■ jq) in Landau gauge within the framework of global color 
model. The result for the dynamical gluon mass is within the range given by other inde- 
pendent determinations. The obtained mixed Landau gauge condensate (qgA-^/q) is clearly 
dependent on the definitions of the condensates. We show that the consistent result may 
be obtained when the same definitions are used. 
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1. Introduction 

The vacuum condensates are believed to play an important role for characterizing the non- 
perturbative aspect of QCD. The well known gauge invariant condensates such as (qq) and 
(a s G 2 ) are proved to be crucial in the application of SVZ sum rules. Recently, there has been 
growing evidence for the existence of the dimension 2 condensate (Af^A®) in SU(N C ) pure 
Yang-Mills theory and QCD. Lavelle and Schaden [Q was responsible for first introducing the 
non vanishing (A^A^) to study the gauge dependent Green functions by making use of the 
operator product expansion (OPE). The recent evidence for the importance of the condensate 
(A^Afj) based on OPE can be found in [21 El II]- The non- vanishing gluon pair condensate gives 
rise to the dynamical gluon mass which is determined by 

m 2 9 = {N 2 N _ C 1)D ( a + 5 " D )(9 2 K A V> C 1 ) 

where a stands for the Lorentz covariant gauge fixing parameter and D is the space-time 
dimension. 

Early implications of the gluon pair condensate have been explored by Celenza and Shakin pj. 
Naively, the operator A 2 is not gauge invariant. The importance of the condensate raises the 
question of gauge invariance and leads to a number of papers which discuss how this condensate 
could play a role in gauge invariant formulation [HJ 13 151 El E2] • Kondo jS] introduced a BRST 
invariant and Anti-BRST invariant combination of dimension 2 operator 



= ^-1 d D xtr 



A tl (x)A ll (x) - iaC(x)C(x) 



where C(x) and C(x) are the Faddeev- Popov ghost and O is the volume of the D-dimensional 
space-time. The vacuum expectation value (VEV) of this dimension 2 operator reduces to the 
gluon condensate in Landau gauge with a = which has a definite physical meaning [7j . 

The role of the dimension 4 mixed quark gluon condensate (gq^-A^X^q) has been studied 
in I12j in the Lorentz covariant gauge and it has been shown that this condensate has the 
contributions to the quark propagator and the gluon propagator. This mixed condensate is 
clearly gauge non-invariant. Recently, the authors of Ref. jl] have evaluated the values of the 
gluon pair condensate and the dimension 4 mixed condensate from the quenched lattice results 
for the quark propagator in the Landau gauge and claimed that it is the first determination of 
the mixed condensate (gq^-A^X^q). This has motivated us to investigate both of these Landau 
gauge condensates from the Global Color Model. 

2. Formalism 

The Global Color Model(GCM) is a quark-gluon filed theory with the action 

S GC M[q,q,A;} = J (qfr-d + m-igA^iJq+^AlDrfiidW) (3) 
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and the generating functional 

Z[J, V, V }= J DqDqDAexp(-S GC M % q, A J] +rjq + QV + ^). (4) 

The essence of GCM is that it models the QCD local gluonic action J F^ V F^ V which has local 
color symmetry, by a highly nonlocal action which has a global color symmetry. To quote 
Cahill and Gunner [T^]: "There is an Infrared(IR) saturation effect which in conjunction with 
the dynamical breaking of chiral symmetry, appears to suppress details of the formal color 
gauge symmetry of QCD." The main aspects of GCM have been reviewed in |13 [ I14 | IT5]. 

In Euclidean metric, with {7^,7^} = 25^ v and 7+ = 7^, the inverse of the dressed quark 
propagator in the chiral limit has the decomposition 

G~ l {p) = ij ■ p + S(p) = ij-pA(p 2 ) + B(p 2 ), (5) 

where stands for the dressing self-energy of quarks. Within the GCM formalism, the 

quark self-energy is determined by the rainbow-ladder truncated quark DSE 

E(P) = I / -0^9 2 D^(p - q)^G(q) lv , (6) 

where g 2 D fMU (p — q) is the effective gluon propagator. In Landau gauge with g 2 D flu (p) = 
{fifjtu — ^jr-)g 2 D{p 2 )i the nontrivial solution to Eq.® which is characterized by B{p 2 ) ^ is 
determined by the two coupled integrable functions 



2 \ -1 \ 2 4 f dA q 2 n / m , n q ■ (p- q)p- (p - q)s Mq 2 '' 



while the trivial solution to Eq.© which is characterized by B{jp 2 ) = is determined by 

( A i(2\ in 2 4 / d4< 7 2 n / \t L0 9- (P-?)P- (P~?)n 1 rn N 
(P ) " 1)P = 3 y ^ - Q)(P ■ 1 + 2 )^^5) • o) 

Within the GCM formalism, the "Nambu-Goldstone" phase is described by the nontrivial quark 
propagator since chiral symmetry is dynamically broken with the no-vanishing constituent 
quark mass M(p 2 ) = B(p 2 )/A(p 2 ) in the chiral limit and the dressed quarks are confined 
since it does not have a Lehmann representation. The "Wigner" phase is described by the 
corresponding trivial quark propagator G'(p) = j^jij^rj with neither chiral symmetry breaking 
nor confinement. The quark propagators G(x — y) and G'[x — y) can be seen as the expectation 
values of the operator Tq(x)q(y) over the physical vacuum state \V) and the perturbative 
vacuum state \P) respectively within GCM formalism. 
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Since the functional integration over the gluon field A® is quadratic according to ©, one 
can perform the integration over gluon field analytically. Taking the quark color current = 
Wq^flnQ as the external source term for the gluon field A a ^ we have the typical gaussian 
integrations 

' VAe~^ AD ~ lA+jA = e^ Dj 



J VAAe~\ AD lA+ i A = (jD) e y^ (10) 

VAA 2 e -\AD^A +j A =[D + ( jD f ]e hjDj_ 
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Through above technique, the gluons vacuum averages are replaced by the quark color current 
together with the effective gluon propagator D. This method is very similar to the deter- 
mination of the VEVs of the QCD operators involving the gluon fields in the instanton liquid 
model, where by integrating over the instanton coordinates, one derives an effective quark ac- 
tion of the form of a Nambu- Jona-Lasinio model ^0] • This method was first used by Meissner 
|17j to calculate the gauge invariant dimension 5 mixed quark-gluon condensate {gqaGq) in 
GCM. In the mean field level, it is straightforward to calculate the vacuum expectation value 
(VEV) of any quark operator of the form 

°n - ^AgUj^A^fe) • • • (^a£UJ, (11) 

where stands for an operator in Dirac, color and flavor space. The expression for the VEV 
of the operator O n takes the form 

{o n ) = (-i) n E(-) p [41^-- A £L^u K1 )---^ PW ]> (i2) 

V 

where p stands for the permutation of the n indexes. According to (|12|) . one can obtain the 
familiar expression for the quark condensate 

m = ~^L ' bA*(s)+ W (13) 
The expression for the four quark condensate (qA^qqA^q) can also be obtained within this 
formalism which is consistent with the result based on the vacuum saturation approximation. 

Note that there is UV divergence in the calculation of the vacuum condensates based on 
Eq. (|12[) . The UV divergence of Eq. (jl2j) can be traced back to the UV behavior of both the 
scalar quark self energy function B(p 2 ) and the vector quark self energy function A(p 2 ). The 
ultraviolet behavior of the constituent quark mass M(p 2 ) = B(p 2 )/A(p 2 ) takes the form in the 
Landau gauge JH] 

M(p 2 ) = 2 -^ , (14) 

3 p 2 ( i M ^ D i- 7ro ' 

^QCD 
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where "y m = 12/(33— 2Nf) is the mass anomalous dimension, with Nf the number of light-quark 
and (qq}° the renormalization-group-invariant condensate. It is easy to see that the Ea. (fT3*|) 
is logarithmic divergent according to the UV behavior of B(p 2 )/ A(p 2 ). In Ref |19) . the quark 
condensate is well-defined and the UV divergence problem of the quark condensate is conquered 
by multiplying a renormalization constant. Within the framework of GCM, Eq. (jl3|) will be free 
of UV divergence if the chosen effective gluon propagator has finite momentum range. Since the 
vacuum condensates are mainly determined by the low and medium momentum properties of 
the Schwinger functions of QCD, ignoring the effects from the hard gluonic radiative corrections 
will have little impact on the determination of these condensates. Note that the definition of the 
vacuum condensate in Ref ^7j; which used a finite cutoff on the integral, is incorrect since the 
obtained values for the quark condensate and the mixed quark-gluon condensate are distinctly 
sensitive to the cutoff. Without the cutoff, exploring of the thermal properties of g(qaGq) 
based on GCM formalism was made in Ref. 20 . 

Even in the case of not taking into account the asymptotic behavior of the effective gluon 
propagator, the UV behavior of the function A(s) can also lead to the UV divergence problem 
in the calculation of the condensates based on Eq. (|12|) . In Ref |21| , the calculation of the mixed 
tensor susceptibility was performed based on (|T2|) where the UV divergence caused by the UV 
behavior of the function A(s) was cancelled by subtracting the corresponding perturbative 
contribution which was evaluated by the trivial quark propagator G'(p). The UV divergences 
induced by the function A(s) also appear in the calculation of the Landau gauge condensates 
in this letter. Actually, this kind of divergence is due to no subtraction of the perturbative 
contribution in the calculation since in this case the Ea.(|12|) is also divergent if one substitutes 
G'{x — y) for G(x — y). Because the condensates reflect the nonperturbative structure of 
QCD vacuum, the perturbative contribution to Eq. (|12|) must be subtracted. In the framework 
of GCM, the perturbative contribution can be self consistently evaluated by replacing the 
nontrivial quark propagator G in Ea.(|12|) with the corresponding trivial quark propagator 
G' . As a matter of fact, the trivial quark propagator has been used extensively to play the 
role of the perturbative dressed quark propagator in the study of thermal properties of QCD 
within the DSE formalism, where the bag constant was defined as the difference of the pressure 
between the true QCD vacuum and the perturbative QCD vacuum, which were evaluated by the 
Nambu-Goldstone solution and the Wigner solution to the quark propagator, respectivelv|22j . 

Therefore, the formula for the calculation of the VEV of operators in terms of quark fields 
should be changed by subtracting the corresponding trivial contribution 

<:£>„:) = (-l) n E(-) ? ^!r- A Sll G Hi f(1 ,'--^ W ] 
p 

-i-ir B-n A 5;L • • • a£L<W • - G UJ- (15) 
p 
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In fact, for the cases of the determination of the quark condensate and gauge invariant mixed 
quark gluon condensate, there are no needs to subtract these trivial terms which are rigorously 
zero due to B(p 2 ) = 0. 

Since the GCM is not renormalizable and the adopted effective gluon propagator in our 
calculation has a finite range in the momentum space, the scale at which a condensate is 
defined in our method is a typical hadronic scale, which is implicitly determined by the chosen 
effective gluon propagator D(s) and the corresponding solution of the quark gap equations. In 
fact, the situation is very similar to the calculation of the condensates based on the instanton 
liquid model where the the scale is set by the inverse instanton size. 

3. Evaluation of Landau gauge condensates 

With the above preparation, the gluon pair condensate and the dimension 4 mixed quark- 
gluon condensate can be evaluated in the framework of GCM. Using Ea. (|lUj) and Ea. (|15|) , we 

get the condensate g 2 {A a (0)^(0)) 



g 2 (A a Jx)A a Jx)) 



(g 2 D™(x-x))- j 'dV f d%'g 2 D%\x-x')g 2 D%,{x-y') 



x(: ^{x')^^—q{x')q{y')^ v ,'-^-q{y') :) 



2 



(16) 



where 



A" X b 
(■ q{x')-l^—q{x')q(y')-i u >—q{y') 



2 



-tr[ 7 ^G'(x' - y >) lv ,^G'(y' - x')]. (17) 

From Eq. 1)16(1 , one can see that the contribution to the gluon pair condensate is divided into two 
parts in the mean field level: the contribution from the effective propagator and the contribution 
from the four quark condensate. After Fourier transformation, we have 

''" '" d 4 p d A q 2 



(2vr) 4 (2tt)' 



: g z D(p - q)g 2 D(p 



x{l2B(p 2 )B(q 2 )Z(p 2 )Z(q 2 ) + 



4p ■ q + 8 



(p - q) ■ p(p - q) 



(p - q) 2 

A(p 2 )A(q 2 )Z(p 2 )Z(q 2 ) - A'{p 2 )A'{q 2 )Z\p 2 )Z\q 2 ) 

oo pi 

dt / 





2W0 dssa2D{s)+ ^L ds 



dxstyl 



-1 



x [g z D(s + t - 2xVst)] z < 3Z(s)Z(t)B(s)B(t) + 



(s — yfsix){^fsix — t) 



s + t — 2^/six 



+Vstx 



Z(s)Z(t)A{s)A(t) - Z'{s)Z'{t)A'{s)A'{t) 



(18) 



6 



with Z{p 2 ) = l/{A 2 {p 2 )p 2 + B 2 {p 2 )) and Z'(p 2 ) = l/(A' 2 (p 2 )p 2 ). 

Applying the same method, we obtain the expression for the mixed condensate in Landau 
gauge 



igiqAi—^q) 



-4 J d*x'g 2 D lxv {x-x'){tr D [ 1 „G{x-x') lv G{x' - x)\ 
-tr D [j^G\x - x')-f v G'{x' - x)]j 



-12 



d 4 p 



B 2 (p 2 



+ 



A{p 2 )(A{p 2 ) - l)p 2 



(2vr) 4 \lA 2 {p 2 )p 2 + B 2 {p 2 ) A 2 (p 2 )p 2 + B 2 {p 2 



A'{p 2 



A'(p 2 



3 
4^2 



dss 



A(s)s 



(A 2 (s)s + B 2 (s)) A'(s) 



Note that the quark DSE 



S(aj -y) = -g 2 D^ u (x - y)7^G(x - y)-y u 



(19) 



(20) 



has been used again in above derivation, which strongly simplifies the evaluation of (j!9j) . From 
(|19j). one can see that the mixed condensate vanishes when QCD undergoes a phase transition 
from the Nambu-Goldstone phase to the Wigner phase with B — > and A — > j4' . 

As mentioned above, the reduced version of Eqs.(|lHjl and ()19|) without the corresponding 
subtraction term will be UV divergent. The UV divergence arises from the behavior of A{p 2 ) 
at large momentum region. In the GCM formalism, A(p 2 ) decreases with increasing p 2 until 
it takes the value A(p 2 ) = 1 at p 2 = oo. This property can be illustrated by a simple analytic 
confining infrared-dominant (ID) model 

g 2 D^{p -q) = {2^5\p - q){5^ - (21) 



4 - • • ■ p 

which was proposed in Ref j22j. The ID model has been used extensively in the literature, 
which can exhibit many of the qualitatively features based on more complicated ansantz. The 
Nambu-Goldstone solution for this model is 

(22) 



A(p 2 ) 
B{p 2 ) 



P 2 < V ~, 



1(1 + a/1 + 



P 2 <\, 



\ yjrj 2 - Ap 2 , 

1 0, otherwise 
The corresponding Wigner solution takes the form 



otherwise, 

2 



1 



sV) = o f A'{p 2 ) = t{i + Ji + 2 4- 



(23) 



(24) 
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We can see that the function A(p 2 ) is identical to A'(p 2 ) in thep 2 > 2L. region and monotonously 
decreases with p 2 until it reaches 1 at oo. Using the ID model, the A(p 2 ) induced UV divergence 
without subtraction can be easily shown. For example, due to A(s) — 1 ~ - for s — > oo according 
to (|22[). the second term of right hand side of Eq. (|19|) is clearly divergent if we do not adopt 
the subtraction procedure. From Eg. 1)221) and Ea.(|24|). it is also shown that the contribution 
of the last term in Eg. (|19|) to the mixed condensate is mainly determined by the low energy 
behavior of the A(p 2 ) and A'(p 2 ). 

The chosen model ansatz for g 2 D(s) in this letter is the Gaussian type model 

g 2 Df u (q) = 4tt 2 a4 exp(-4)5 a V(</) (25) 

p 0J Z 0J Z 

where the parameter A sets the strength of the interaction and u is the range parameter at 
which the scalar function in the parameterization is maximal. This Gaussian type form was 
supposed to represent a sensible hadron model in |24j and the parameter oj was envisaged 
to have a value of several hundred MeV which sets the interaction scale. In our work, these 
parameters are constrained by giving the reasonable constituent quark mass M q (q 2 = 0) and 
the pion decay constant fixed at 87 MeV which is more appropriate in the chiral limit rather 
than the Pion's mass-shell value of 93 MeV. To check the sensitivity of our results on the model 
effective gluon propagator, various sets of parameters A and oj were used in our calculations. 

Table 1. Calculated Landau gauge condensates from the Gaussian model in the chiral 
limit with three sets of parameters. The quark condensate and the constituent quark 



mass 


from this model are 


also listed. The units for A, oj, M g , 


M q and f n are 


GeV. 


A 


oj (qq)(GeV 3 ) 


g 2 (AA)(GeV 2 ) 


g{qA- iq )(GeV 4 ) 


Mg M q (0) 


U 


30.0 


0.40 -0.215 3 


9.5 


-1.2* 10" 3 


0.94 0.381 


0.087 


19.0 


0.45 -0.223 3 


9.2 


-1.3* 10" 3 


0.93 0.342 


0.087 


13.2 


0.50 -0.232 3 


9.9 


-1.4* 10" 3 


0.96 0.311 


0.087 



Table 1 shows that the dynamical gluon mass obtained in the Gaussian model is around 
lGeV. This value is compatible with the estimates obtained by other approaches, which are 
within the range 0.5 — 1.5GeV |25j . Actually, from our numerical studies, the main contribution 
to the gluon pair condensate comes from the first term of right hand side of Eq. ()16|) . This is 
shown in Table 2 with the dynamical gluon masses M g \ which is determined only from the 
first term of the Eq. (|16[) . Comparing M g \ with M g , we get the conclusion that the gluon pair 
condensate is mainly determined by the low energy part of g 2 D(p) and the correction from 
the four quark condensates can be neglected. This is reasonable since in GCM the g 2 D(p) is 
interpreted as the effective gluon propagator. Our results further confirms that the Gaussian 
model is surely an effective quark-quark interaction form in the low momentum region. In 
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addition, this also provides a constraint to the form of g 2 D(p) in GCM, which should give a 
reasonable dynamical gluon mass. 

Neglecting the four quark condensate contribution, the gluon pair condensate from the ID 
model takes the simple form 

g 2 (A 2 ) = 2N cV 2 . (26) 

In the case rj /2 = m 2 p = 0.59GeV 2 |23], the gluon pair condensate takes the value 7.1GeV 2 and 
the corresponding gluon mass is 0.82GeV which is also in agreement with the estimates from 
other methods. 

In Table 1, the magnitude of the obtained dimension 4 mixed condensate is much smaller 
in comparison with the recent results based on the quenched lattice results for the quark 
propagator in Landau gauge 4 . Within the ID model with rj 2 /2 = 0.59GeV 2 , the mixed 
condensate takes the value — 0.0003GeV 4 . However, the magnitude of the quark condensate 
— (qq) 1 ^ obtained in ID model with the same parameter is 118MeV which is also smaller 
in comparison with the standard value 250MeV. If we simply lift the the value of r/, for 
example, taking r/ = 1.8GeV, we get that g{qA ■ 79) = -0.0018GeV 4 with -{qq) 1/3 = 184MeV. 
Therefore, the mixed condensates obtained in both the ID and Gaussian models are at the 
same order of magnitude. According to [TT 1121 Hj. the coefficient of 1/Q 4 in the expansion 
of the vector quark self-energy function at sufficiently large values of Q is determined by the 
combined dimension 4 condensates 

a s (qgA. iq )-^(^G 2 ), (27) 

2 ( TT 

where the gluon condensate takes the value (^-G 2 ) ~ O.OlGeV 4 . Multiplying a s (fi) ~ 0.5 to 
the mixed condensate where [i stands for a typical hadronic mass, we find that a s {qgA ■ 'jq) 
is around — 0.0007GeV 4 in our calculation. In comparison with the magnitude of (^-G 2 ), we 
find that the main contribution to Eq.(J57J) is determined by the gluon condensate. This result 
is in contrary to the one obtained in Ref 4 . The obtained value for a s (qgA ■ 75) in Ref [1] 
is (—0.11 ± 0.03)GeV 4 , comparing to which the gluon condensate contribution to (|27|) can be 
neglected. 

Table 2. The comparison of the contributions to the gluon pair condensate from 



the effective gluon propagator and the four quark condensate in Eg. 1)16(1 . 



A(GeV) 


w(GeV) 


g 2 {AA) 1 (GeY 2 ) 


g 2 (AA) 2 (GeV 2 ) 


M sl (GeV) 


30.0 


0.40 


9.22 


0.28 


0.93 


19.0 


0.45 


9.35 


-0.13 


0.94 


13.2 


0.50 


9.90 


-0.05 


0.96 



According to the quark motion equation (7^^ + m — i-^-A^^q = , we expect that the 
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condensate (q^^d^ 
reduced to 



m 



■^-A'^'y^q) should take zero. In the chiral limit, this relation is 



ya 



(28) 



Actually, the evaluation of the Landau gauge condensate {q^^d^iq) is very simple in GCM. 
Repeating the procedure above, including the subtraction, we get 

-(?WtA<|(x)) = tf{tr[G(x'-x)^}-trlG'(x'-x)^}}\ x ,= x 



12 

3 
4^2 



d 4 p 



A(p 2 )p 2 



dss 



A 2 {p 2 )p 2 + B 2 {p 2 ) 
A(s)s 



A'(p 2 



A 2 (s)s + B 2 (s) A'(s) 



(29) 



It can been see that this is just the Ea. (|19|) . However, without the subtraction, Eq. (|19|) reduces 
to 

B 2 (s) A(s)(A(s) 



3 



dss 



+ 



Is 



A 2 (s)s + B 2 (s) A 2 (s)s + B 2 (s) 
which is negative and divergent since A(s) > 1 for finite s, while Eq, (|29j) reduces to 



3 
4^2 



dss 



A(s)s 



A 2 (s)s + B 2 {s) 



(30) 



(31) 



which is positive and divergent. Therefore, there exists a large difference between these two 
expressions even in sign if we do not adopt the subtraction mechanism. This shows that it 
is reasonable to introduce the subtraction mechanism in our definition. For the massive case, 
according to the PCAC relation, the value of m(qq) is around — 0.00008GeV 4 , which suggests 
that the Landau gauge dimension 4 mixed condensate is an order of magnitude larger than the 
m(qq) condensate. 

Here we still want to stress that the discrepancy of the estimates of the mixed quark gluon 
condensate arises from the different definitions of condensates. In Ref4_ and those related to it, 
the condensates are defined as the coefficients of expansion terms in A(s). Actually, by taking 
the same definition, similar results can also be obtained within Global Color Model. Consider 
the ID model described by Eqs. (|22j) and (|2T ~ 



A(s) 



1 



l + y/T+ 2r] 2 /s 



large-s ffl 

2 s 



-rj 2 ^ 1 
_~2~ 



+ 



With the conventions of Ref . [I] one would read 



a s {A 2 



Ncrf 
vr 2 



a s {qgA ■ jq) 



9N r " 7T 



(-G) 



4 

3^ 



(32) 

(33) 
(34) 
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Note that Eq.® is identical to Eq.©. Substituting rj 2 /2 = 0.59GeV 2 into Eqs.jHHJ) and ® 

one gets 

a s (A 2 ) = 0.57 GeV 2 , (35) 
a s (qgA. 1 q)--^-(^G 2 ) = -0.15 GeV 2 , (36) 

9N C 7T 

which are in agreement with the results of Ref . 4 . 

The small value of the mixed quark gluon condensate obtained based on the definition (|12|) 
can be tracked back to subtracting the coefficients 1/s and 1/s 2 in the expansion of A'(s) 
evaluated in the Wigner phase from those obtained with A(s) calculated in the Nambu phase. 
In the case of sufficiently large value of s 

A{s) = 1 + c n\ + c n1_ A {s) = i + C w l _ + C w±_ (37) 
s s z s s z 

one gets 

ocfaA ■ iq ) - ^-(^G 2 ) = 4 - c™. (38) 

The fact that the vector self-energy function of the quark propagator is not too sensitive to the 
phase in large momentum, 1.6., C<2 ^ C<2^ i Cclll qualitatively explain why the l.h.s of Eq. (j38|) is 
small according to our definition. 

4. Discussion and conclusion 

We have evaluated the gluon pair condensate and the dimension 4 mixed quark gluon con- 
densate in Landau Gauge within the framework of GCM. To avoid divergence and subtract 
the perturbative contribution, two measures are taken in our calculations. The input effective 
gluon propagator adopted in our calculation has finite range in momentum space, which can 
suppress the divergence induced by the UV behavior of the scalar quark self-energy function. 
This is reasonable since the condensates are mainly determined by the low energy properties 
of QCD. To subtract the perturbative contribution, the trivial solution or the Wigner solution 
to the quark DSE was used self-consistently in our formalism. 

With above method, the contributions to the gluon pair condensate was decomposed into 
two parts: the contribution from the effective gluon propagator and the contribution from the 
four quark condensate. From our numerical study, we find that the correction to the gluon pair 
condensate from the four quark condensate is very small and can be neglected. The obtained 
gluon pair condensate and the corresponding dynamical gluon mass are compatible with other 
estimations in the literatures. The expression for the dimension 4 mixed quark gluon condensate 
is very simple in our formalism. This expression has been verified by using the quark motion 
equation. The value of the mixed condensate (gqA ■ ^q) is an order of magnitude larger than 
the condensate m(qq). 
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We also explicated that the discrepancy of the estimates of the Landau mixed quark gluon 
condensate between our formalism and others arises from the different definitions of the con- 
densates. If we take the same definition as in Ref.Jl], the obtained values of both the gluon 
pair condensate and the mixed condensate in GCM are consistent with the corresponding ones 
extracted from Lattice simulation results. That means when the same definitions are used, the 
GCM/DSE formalism is completely consistent with the analysis based on the lattice results. 
Note that due to the ambiguity in the definitions of the condensates, there exists no natural 
or unique definitions but only clear definitions; when the same definitions are used, all well- 
constrained models should yield the same result. This fact is made very clear for the quark 
condensate (qq) |2*o] . 
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